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Abstract. The T c k location of Pre-Main Sequence (PMS) evolutionary tracks depends on the treatment of over-adiabaticity 
(D'Antona & Mazzitelli 1994, 1998). Since the convection penetrates into the stellar atmosphere, also the treatment of convec- 
tion in the modeling of stellar atmospheres will affect the location of the Hayashi tracks. In this paper we present new non-grey 
PMS tracks for T e s > 4000 K. We compute several grids of evolutionary tracks varying: i) the treatment of convection: ei- 
ther the Mixing Length Theory (MLT) or Canuto et al. (1996, CGM) formulation of a Full Spectrum of Turbulence; ii) the 
atmospheric boundary conditions: we use the new Vienna grids of ATLAS9 atmospheres (Heiter et al. 2002a), which were 
computed using either MLT (with a = A/H p = 0.5) or CGM treatments. For comparison, we compute as well grids of models 
with the NextGen (Allard & Hauschildt 1997, AH97) atmosphere models, and a 1 M Q grey MLT evolutionary track using the 
a calibration based on 2D-hydrodynamical models (Ludwig et al. 1999). These different grids of models allow us to analyze 
the effects of convection modeling on the non-grey PMS evolutionary tracks. We disentangle the effect on a self-consistent 
treatment of convection in the atmosphere of the wavelength dependent opacity, from the role of the convection model itself in 
the atmosphere and in the interior. We conclude that: i) In spite of the solar calibration, if MLT convection is adopted a large 
uncertainty results in the shape and location of PMS tracks, and the MLT calibration loses sense, ii) As long as the model of 
convection is not the same in the interior and in the atmosphere, the optical depth at which we take the boundary conditions 
is an additional parameter of the models. Hi) Furthermore, very different sub-atmospheric structures are obtained (for MS and 
PMS stellar models) depending not only on the treatment of convection, but also on the optical depth at which the boundary 
conditions are taken, iv) The comparison between NextGen based models and ATLAS9 based models shows that in the T sff 
domain they have in common (4000-10000K) the improved opacities in NextGen atmosphere models have no relevant role on 
the PMS location, this one being determined mainly by the treatment of the over-adiabatic convection, v) The comparison be- 
tween theoretical models and observational data in very young binary systems indicates that, for both treatments of convection 
(MLT and CGM) and for any of the atmosphere grids (including those based on the 2D-hydrodynamical atmosphere models), 
the same assumption for convection cannot be used in PMS and MS: either the models fit the MS - and the Sun in particular- 
or they fit the PMS. Convection in the PMS phase appears to be less efficient than what is necessary in order to fit the Sun. 
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1. Introduction 

The experimental data presently available for young objects 
in star-formation regions need to be compared with theoreti- 
cal evolutionary tracks in order to be correctly interpreted in 
terms of age, mass and chemical composition. As discussed 
in D'Antona (2000) and the references therein, the HR dia- 
gram location of a star during its Pre-Main-Sequence (PMS) 
evolution is very sensitive to physical inputs such as low- 
temperature opacity, equation of state, rotation, atmosphere 
model, and convection treatment. 

During the last years, a lot of work has been done to im- 
prove the knowledge of low-temperature opacities and to in- 
clude them in the modeling of stellar atmospheres. The variety 
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of theoretical evolutionary tracks available today (new convec- 
tion treatment with grey boundary conditions; classic convec- 
tion with non-grey atmosphere models...) has contributed to 
create some confusion about the effect of the different physi- 
cal inputs on the results. Our aim is to extricate the different 
roles of non-grey atmospheres and of convection in evolution- 
ary track computations. As shown by Montalban et al. (2001), 
that can be done only using models with a self-consistent treat- 
ment of convection both in the atmosphere and in the interior. 

In stars on the cool side of the HR diagram, convection is 
deep and defines the "envelope" portion of the star in which 
the dominant mode of energy transport is convection. In par- 
ticular, Pre-Main Sequence (PMS) models are largely convec- 
tive and — due to the low gravity and low temperature — the 
convection can be over-adiabatic in extended regions of the 
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stars. Consequently, the shape of the corresponding evolution- 
ary tracks depends on the efficiency of the convective transport 
of energy. 

The convective zone increases with decreasing r e fj and it 
is gradually shifted towards larger depths, so that the structure 
of the surface layers is less and less affected by the convection 
as T e g decreases. For temperatures lower than 4700K, however, 
the convection zone rises again due to the dissociation of H2. 
On the other hand, for given T e ff and metallicity, and decreas- 
ing gravity, the convective flux decreases because of the lower 
density and, therefore, the over-adiabatic region in the atmo- 
sphere is more extended. Consequently, in cool stars, a large 
part of the photosphere actually forms the uppermost portion 
of the convective envelope. 

Several facts indicate that the classical grey atmosphere 
approximation adopted in stellar computation is not valid for 
cool stars: the diffusion approximation is not valid up to t = 
2/3 (Morel et al. 1994), and the convective transport of en- 
ergy has to be taken into account in the atmosphere as well. 
Furthermore, the introduction of the frequency dependence of 
opacity modifies the onset of convection, and the grey approx- 
imation produces errors in the theoretical T t g and in the col- 
ors of low mass stars (Baraffe et al. 1995; Baraffe & Chabrier 
1997). Therefore, in the computation of stellar models, the 
boundary conditions (BCs) at the surface must be provided by 
non-grey modeling of stellar atmospheres, taking into account 
the dependence of opacity on frequency and an adequate treat- 
ment of convection. 

Unfortunately, for the time being, only very simple local 
convection models are available for routine computation of ex- 
tended grids of model atmospheres, while detailed numerical 
simulations are still unaffordable for applications such as stel- 
lar evolution that requires the calculation of many thousands of 
individual model atmospheres over the HR diagram. 

The "standard model" of convection in stellar evolution is 
the mixing length theory (MLT, Bohm-Vitense 1958), where 
turbulence is described by a relatively simple model that con- 
tains essentially one adjustable parameter, the mixing length: 
A = aH p (Hp being the local pressure scale height and a an 
unconstrained parameter). Canuto & Mazzitelli (1991, CM) 
and Canuto et al. (1996, CGM) made available an alternative 
model (Full Spectrum Turbulence -FST- models), that over- 
comes some of the problems of MLT, while keeping a low 
computational cost. The main characteristics and differences 
of these models are described in Sect 2. 

D'Antona & Mazzitelli (1994, DM94) published four sets 
of PMS evolutionary tracks that allowed to study the effect of 
low-temperature opacity and treatment of convection on stel- 
lar evolution. These tracks were computed with both kinds of 
convection treatments: MLT calibrated on the Sun, and the FST 
model with the CM formalism. The results showed that the T e g 
location of Hayashi tracks has a strong dependence on the treat- 
ment of convection. DM94's models, widely used and tested 
with observations, were revised in D'Antona & Mazzitelli 
(1997, 1998 — DM97, 98) by introducing several improve- 
ments in the micro-physics (updated opacity tables and equa- 
tion of state) and in the macro-physics (the updated FST treat- 
ment of convection by Canuto et al. 1996). The DM97 models 



were still employing grey atmospheric BCs, but the problem 
of matching convection in the interior and in the atmosphere 
occurs already in these models. In fact the DM97,98 models 
make only an exploratory approximation for defining the con- 
vective scale length if convection penetrates in the atmosphere, 
namely, they do not include the atmospheric convective depth 
in the computation of the scale. With this choice, convection in 
the sub-photospheric layers is less efficient, and very low mass 
stars are cooler by up to ~ 150 K. 1 DM98 suggested that since 
the over-adiabaticity is present in low gravity / low tempera- 
ture atmospheres, the FST convection treatment would modify 
convection also in the atmosphere, therefore it would be very 
important to include a revised treatment of convection in model 
atmospheres. 

Recently, Heiter et al. (2002a) have published new atmo- 
sphere model grids based on Kurucz's ATLAS9 code (1993). 
They performed calculations for different treatments of non- 
adiabatic convection: MLT (a = 0.5), CM, and CGM (for a 
smaller range of parameters, MLT — a = 1 .25 — was used as 
well). These choices for modeling convection are extensively 
motivated in Sect. 2 of that paper. As one of their main results, 
the authors conclude that MLT models with a small mixing 
length parameter (e.g., a ~ 0.5) and FST models are equivalent 
in the atmospheric region where the observed flux originates. 
Both treatments predict a low convective efficiency for these 
layers. The deep atmospheric structures, however, are differ- 
ent, and each T(r) relation represents stars which differ in ra- 
dius and luminosity, hence the PMS tracks will also be differ- 
ent. Furthermore, if MLT is adopted, and a low a-value is used 
in the atmosphere (a- dtm ) (a atm =0.5, as adopted in Heiter et al. 
2002a, or ar a tm=l as in Hauschildt et al. 1999a 2 ), we must com- 
pensate, in order to fit the Sun, for the high over-adiabaticity in 
the atmosphere, by using an a-value in the interior (aj nt much 
larger than a atm ). FST, by construction, reproduces this behav- 
ior: it is very inefficient at the outer boundary of the convec- 
tion zone, and very efficient in the inner layers. As a conse- 
quence, FST has the advantage of simultaneously (and thus 
consistently) fitting, without arbitrary tuning of the parameter 
set (a- dtm , T P h, ai nt ), the Balmer line profiles (Heiter et al. 2002a) 
and the solar radius (Heiter et al. 2002b and this paper). 

In order to analyze the impact of sub-photospheric convec- 
tion on the HR location of the Hayashi tracks, we computed 
solar composition stellar models for masses from 0.6 to 2.0 M 
with several combinations of convection model and boundary 
conditions (Table 1). The physical inputs of these models are 
described in Sect. 3. 

An additional problem to be considered when computing 
stellar models with non-grey boundary conditions is the choice 
of T p h. From this point, the integration of the atmospheric strat- 
ification supplies the values of T e g and log g (or luminosity and 



1 The procedure is extensively explained in DM98. Although this is 
the only case in which FST convection results in a lower global effi- 
ciency than MLT, this result has been ascribed in Baraffe et al. (2002) 
to the general behavior of FST large temperature gradients in the at- 
mosphere. 

2 NextGent atmosphere models named here AH97 were partially 
published in Hauschild et al. 1999a. 
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radius) for a given mass. If the atmospheric integration is con- 
sistent with the physics in the interior, and the diffusion approx- 
imation holds below r p h, the model location in the HR diagram 
should not depend on the choice of r p h. Since grids of atmo- 
sphere models and interiors are usually computed by different 
teams with different aims, quite often the physics (convection 
model, opacity, equation of state, etc.) used in both regions is 
not exactly the same. In Sect. 4 we analyze how these differ- 
ences affect the evolution and structure of 1 M Q star. 

In Sect. 5 we present our new non-grey PMS evolutionary 
tracks for T e g > 4000 K. We have also computed two grids 
of "complete" FST models with metallicity larger and smaller 
than the solar one by a factor of 2. The effect of different chem- 
ical compositions on HR diagram location is shown in Sect. 6. 

In Sect. 7 we compare our new FST Hayashi tracks with the 
available data for young binary stars in star-formation regions. 
In fact, PMS binaries having experimentally known masses are 
the best test for the PMS models, since they should fit both 
masses and display the same age. Finally, the main conclusions 
of this paper are summarized in Sect. 8. 

2. About convection in stellar modeling 

The system of equations describing stellar convection from 
first principles is both non-linear and non-local and its an- 
alytical or numerical solution for the general stellar case is 
not yet available. Even numerical simulations not resolving 
all the spatial scales of the flow are unaffordable for stellar 
evolution modeling, at least for the intermediate future, be- 
cause of the excessive thermal relaxation time scales of stel- 
lar structure modeling. Non-local models such as those pro- 
posed by Xiong (1985), Canuto (1992), and Grossman (1996) 
are a promising and much more affordable alternative with the 
present computational resources. But for the time being, in 
stellar modeling we still stick to local models: the "standard" 
Mixing Length Theory (Bohm-Vitense 1958), and the "Full 
Spectrum of Turbulence" formalisms (Canuto & Mazzitelli, 
1991; Canuto et al. 1996) are two examples. 

MLT treats the heat transport with a one-eddy incompress- 
ible model with compressibility partially accounted for through 
the definition of a mixing length A = aH p , where H p is the lo- 
cal pressure scale height and a is a free parameter. Mimicking 
the spectral distribution of eddies by one "average" eddy (re- 
liable only for high viscosity fluids) has critical consequences 
on the computation of MLT fluxes: Canuto (1996) showed that 
in the limit of highly efficient convection (S » 1, where S is 
the convective efficiency, see e.g. CM for details) MLT under- 
estimates the convective flux (a fact also confirmed by labora- 
tory data, Castaing et al. 1989), and in the low efficiency limit 
(S « 1) MLT overestimates the convective flux. 

Computations show that inside the stars with deep convec- 
tive regions S <K 1 at the top of the convective envelope, just 
below the surface, while at the peak of the super-adiabaticity 
S » 1. In inefficient convection, the convective temperature 
gradient sticks to the radiative one and begins detaching from 
it only when convection becomes efficient. 

Historically, the choice of a mixing length A = aH p is 
in part a consequence of this underestimation of fluxes in 



MLT. Thus, the scale length A was originally (Bohm & Stiickl 
1967) chosen coincident with the distance z from the convec- 
tive boundary, consistently with the Von Karman law for in- 
compressible fluids (Prandtl 1925), but due to the underestima- 
tion of convective flux, models with A = z were not able to 
fit the Sun. 3 Then, A = aH p grew up as standard choice, with 
the free parameter a calibrated by comparing the solar models 
with the actual Sun. Depending on the physical inputs the value 
of a can vary in between 1.5 and 2.2. If we use this prescrip- 
tion for the mixing length, at the peak of the super-adiabaticity, 
aH p » z, and the fluxes, which are oc A 2 , are artificially in- 
creased. 

FST (CM and CGM) models attempt to overcome the one- 
eddy approximation by using a turbulence model to compute 
the full spectrum of a turbulent convective flow. The treat- 
ment of turbulence in CM and CGM is slightly different. 
Nevertheless, in both cases the derived convective fluxes are 
quite distinct from the MLT ones. CM and CGM convective 
fluxes are ~ 10 times larger than the MLT ones for the S » 1 
limit while, for the low efficiency limit (S <k 1), the CM flux 
is ~ 1/10 of the MLT flux, and the CGM one ~ 0.3 times the 
MLT flux. This behavior yields, in the over-adiabatic region at 
the top of a convection zone, steeper temperature gradients for 
FST than for solar-tuned MLT (and, in the FST framework, 
gradients steeper for CM than for CGM). 

Like MLT, FST models assume the Boussinesq approxima- 
tion. Following the physical argument that the Boussinesq ap- 
proximation leaves no natural unit of length other than the dis- 
tance to a boundary (z) and the vertical stacking of the eddies 
due to density stratification, CM and CGM adopt respectively 
Acm = Z, and Acgm = z + a*H p + op . A The second term in Acgm 
is a fine tuning parameter that allows small adjustments, if ex- 
act stellar radii are needed, e.g., in helioseismology. Canuto et 
al. (1996) stress however, that the role of a* in FST models 
is radically different from that of a in the MLT model. FST 
tuning affects, in fact, only layers close to the boundaries since, 
for inner layers, z quickly grows and becomes much larger than 
a*H Pttop . On the other hand, if we wish to attribute a physical 
meaning to a*, the second term in Acgm can be interpreted 
as an "overshooting term" representing the observed fact that 
convection penetrates in the stable region, and thus the scale 
length cannot decay to zero right at the depth where V ra d = V a d 
(Schwarzschild criterion of stability). In any case, a* in these 
models is not properly used as an overshooting parameter but 
only to define the exact scale length. 

The FST fluxes have also been used combined with other 
definitions of the mixing length. Bernkopf (1998), for instance, 
uses the CM version of the FST fluxes combined with a mixing 
length A = aH p and with a < 1 . 

3 We will use the expression "fit the Sun" throughout this paper 
to indicate that a solar mass track of solar metallicity Z reaches the 
present solar luminosity and radius (and thus the present solar r e ff) 
at the solar age. The fit of the solar luminosity determines the helium 
abundance Y - or better the ratio Z/Y (e.g. Turck-Chieze et al. 1988), 
while the fit of the radius mainly depends on the tuning of the convec- 
tive model. 

4 Note that we adopt a slightly different definition of A for the 
present computation of the interior - see Sect. 3.3. 
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Fig. 1. T(t) relations computed using ATLAS9, a FST treatment of convection with a* - 0.09 (Heiter et al. 2002a) and integrated 
using 288 layers (solid lines), or a MLT one with a = 0.5 and 72 layers (dotted lines and full circles). The lower curves correspond 
to J eff =4000 K and the upper curves to r eff =10000 K, with a step Ar eff =200 K. 



3. The stellar models 

The stellar models presented here were computed using the 
stellar evolution code ATON2.0 (Ventura et al. 1998a). Specific 
details concerning the computation of the pre-main sequence 
and deuterium burning phases are described in Mazzitelli & 
Moretti (1980). We computed models for three metallicities 
[M/H] = -0.3,0.0, and +0.3 with an adopted helium mass 
fraction Y = 0.28. Three different sets of boundary condi- 
tions were considered for solar composition: two of the new 
grids of ATLAS9 atmospheres from Heiter et al. (2002a), MLT 
(a = 0.5) and CGM grids, and, for comparison, the non-grey 
atmosphere models by AH97. In the first case the considered 
masses span the interval 0.6 to 2.0 M , in the second one that 
from 0.6 to 1.5 M . For the other two chemical compositions 
only CGM models were computed. In addition, we computed 
a grey MLT evolutionary track for 1 M using for a the val- 
ues provided by the calibration from 2D-hydrodynamical at- 
mosphere models (Ludwig et al. 1999). 

3.1. Equation of state and opacities 

A complete description of the equation of state (EOS) of our 
code is given in Montalban et al. (2000). For the present mod- 



els, the thermodynamics is from Rogers et al. (1996), for five 
different hydrogen abundances. At temperatures T > 6000 K 
we adopt the OPAL radiative opacities (F) (Roger & Iglesias 
1996, for the solar Z-distribution from Grevesse & Noels 
1993). In the high-density (p) regions the opacities are linearly 
extrapolated (log k vs. log p), and harmonically added to con- 
ductive opacities by Itoh & Kohyama (1993). At lower temper- 
atures we use the Alexander & Ferguson's (1994) molecular 
opacities (plus electron conduction in full ionization) for the 
same H/He ratios as in the OPAL case. 

3.2. Atmospheric structure and boundary conditions 

The new grids of ATLAS9 atmospheres by Heiter et al. (2002a) 
introduced several improvements in comparison with previous 
model grids published by Kurucz (1993, 1998) and Castelli 
et al. (1997). Most of them are related to a finer grid spac- 
ing (Ar e ff , A log g, but also vertical resolution) which allows 
a more accurate interpolation within the grids. For the conve- 
nience of the reader we recall the main characteristics of these 
grids: 

- Three different models of convection: MLT(a = 0.5), FST 
according to CM, and FST CGM. 
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- For FST models, an increase of the vertical resolution in 
the atmospheric integration from 72 to 288 layers ranging 
from logT Ross = -6.875 to logT Ross = 2.094 (where t Ross 
is the average Rosseland optical depth). 

- Effective temperature range: 4000-10000 K, with Ar eff = 
200 K. 

- Gravity range: logg from 2.0 to 5.0, with Alogg = 0.2. 

- Chemical composition: sets for [M/H]=-2.0, -1.5, -1.0, - 
0.5, -0.3, -0.2, -0.1, 0.0, +0.1, +0.2, +0.3, +0.5, +1.0 are 
available. 

From the atmosphere models we built the tables of 
BCs corresponding to a fixed r p h. These tables contain, for 
each (r e ff,log g,[M/H]), the following quantities: temperature 
r(r p h), pressure P(T p h), geometrical depth z(r p h), geometrical 
depth z(to) at the optical depth tq for which V ra( j = V a( j, and fi- 
nally the pressure scale height H ? (tq) = H p top . We built tables 
corresponding to r p h=l, 3, 10, and 100 in order to analyze the 
dependence of the results on r P h. We note, however, that Heiter 
et al. (2002a) advise to switch between model atmosphere and 
stellar envelope at T p h=10, to avoid the discrepancies due to 
the turbulent pressure not included in the atmosphere modeling 
and to reduce the effect of different opacity tables and equa- 
tion of state used at both sides of r p h. The boundary conditions 
for the internal structure are determined by spline interpolation 
of these tables. From the initial T e s and log L we determine P 
and T at the last point of the internal structure (r p h) and the 
derivative of P and T with respect to luminosity and radius. 
An iterative procedure is performed until the P(r p h) and 7\T p h) 
values derived at the boundary r p h from the interior or from the 
atmosphere models converge. In order to take into account the 
effect of the penetration of convection into the atmosphere, the 
thickness of the convection zone in the atmosphere is included 
in the definition of the scale length A in the FST models (see 
the discussion of Eq. below). 

Some parts of the tracks for 0.6 and 0.7 M Q reach T e s < 
4000 K, falling outside the grid of ATLAS9 models. We com- 
pleted the evolution by assuming as boundary conditions the 
values of P and T for the 4000 K grid points. The track is then 
artificial until T e s becomes again larger than 4000 K, but this 
does not introduce large errors in the determination of the fol- 
lowing evolution to the MS. We recall here that the validity of 
the plane parallel approximation for log g > 2, and for the T e g 
domain of interest for our work, has already been shown by 
Hauschildt et al. (1999b) and Baraffe et al. (2002) for the case 
of NextGen model atmospheres. We can hence safely to assume 
its validity for the ATLAS9 model atmospheres as well. 

Concerning the NextGen model atmospheres (AH97), their 
opacities include the contribution of many more molecular 
lines than the ATLAS9 ones, consequently, they are more ad- 
equate to study the T e s region where these opacities dominate 
(roughly, below 4000 K). In AH97 models, the convection is 
treated with MLT and a — 1 . For solar metallicity, the avail- 
able models have T e g from 3000 to 10000 K, and a surface 
gravity from \ogg = 3.5 to 6.0, with Alogg = 0.5; the vertical 
resolution in the atmospheric integration is of 50 layers ranging 
from r = 10~ 7 to 100 (only 9 layers between r = 1 and 100). 




-10 12 
log T 

Fig. 2. Upper panel: over-adiabaticity as a function of the opti- 
cal depth in the atmosphere for two pairs of (r e ff,logg) corre- 
sponding to typical values for a 1 M star during MS and PMS. 
Dashed lines correspond to MLT ATLAS9 (o-atm = 0.5) atmo- 
sphere models, and solid lines to CGM FST ones. Lower panel: 
Entropy jump for the same atmospheres as in upper panel. 

3.2.1 . The treatment of convection in the atmosphere 

The original ATLAS9 code treats the convection by using MLT 
with a - 1.25 (cf. Castelli et al. 1997). In addition, the code 
has the possibility of including a sort of "overshooting". The 
value a=1.25 and the option of approximate overshooting were 
originally adopted to fit the intensity spectrum at the center of 
the solar disk and the solar irradiance. However, Castelli et al. 
(1997) showed that these quantities are much more sensitive 
to the overshooting-on mode than to the value of a itself (see 
also Heiter et al. 2002a) and that, if this treatment is included, 
the H a and H^ profiles of the solar spectrum cannot be simul- 
taneously matched. Several papers on: i) the T e ff determination 
from H a and H^ lines (e.g., Fuhrmann et al. 1993, 1994; Van't 
Veer-Menneret & Megessier 1996; Van't Veer-Meneret et al. 
1998); ii) the theoretical predictions of H a and H^g (Gardiner et 
al. 1999 from ID models, and Steffen & Ludwig 1997 from 
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2D numerical simulations), of (b - y) and c Stromgren in- 
dices (Smalley & Kupka 1999) and of Geneva indices (Schmidt 
1999), and Hi) abundance determinations (Heiter et al. 1998), 
indicate that, even if a ID-homogeneous model cannot explain 
all the spectroscopic and photometric observations, atmosphere 
models predicting temperature gradients closer to the radiative 
one (i.e., those models in which convection is less efficient than 
predicted by MLT-based models with a > 1) are in better over- 
all agreement with observations. As pointed out, however, in 
the above-mentioned papers, not all the color indices over the 
whole r e fl domain can be reproduced with such models either. 
Whereas in other cases, the quality of a fit (such as line profiles 
of H a ) is even independent of any sensible value of a (although 
not of other modifications to the convection treatment such as 
"approximate overshooting"). 

In the new ATLAS9-MLT grids, the convection is described 
as in Castelli (1996) and Castelli et al. (1997) (that is V/A = 
A/6 and y = 0.5), but with a = 0.5 and without "overshoot- 
ing". In the CGM model atmospheres, the convective flux is 
computed as in Canuto et al. (1996), but the characteristic scale 
length is defined as: A = min(z top + a*// p , top ,z b ot + a*# P ,bot) 
where the index "top" and "bot" refers to top and bottom of the 
convective region. For most model atmospheres the difference 
with respect to the original prescription in Canuto et al. (1996) 
is either null or negligibly small, because the temperature gra- 
dient for convection zones which are entirely contained within 
the atmosphere is practically the radiative one, while, for con- 
vection zones extending below the atmosphere, the evaluation 
of A in a pure model atmosphere has necessarily to occur near 
the top of the convection zone (Heiter et al. 2002a). 

Since the effect of varying the fine tuning parameter a* is 
felt only in a small region close to the boundary of the convec- 
tion zone, changing a* - say - by a factor two yields no dif- 
ference in the solar spectrum features, so that its value cannot 
be calibrated from atmosphere models. Therefore, the choice 
of a value of a* - 0.09 was based on the solar calibration per- 
formed by Canuto et al. (1996), but it has to be pointed out that 
this value was determined using grey BCs (Henyey et al. 1965). 

As shown in Heiter et al. (2002a) the temperature struc- 
ture in the deep atmosphere is strongly dependent on the con- 
vection model (Fig. Because of the features of convective 
fluxes described in Sect. 2, the convection modeled with MLT 
(a = 0.5) is very inefficient in the whole atmosphere, while 
the FST models, which are similarly inefficient at the outer 
boundary, rapidly begin to be very efficient in the deeper lay- 
ers of solar-like convection zone. In Fig. |2]we plot the over- 
adiabaticity and the jump of specific entropy corresponding to 
either MLT or FST atmosphere models with two different val- 
ues of r e ff and log g. The stellar spectra are not sensitive to the 
deep atmosphere gradients, and from the point of view of at- 
mosphere models, as long as one uses a low value of a, MLT 
and FST are essentially equivalent. The BCs provided by both 
atmosphere models can become, however, very different as the 
optical depth where the BCs are taken increases. 



3.2.2. Resolution 

Schmidt (1999) showed that a scale length prescription as used 
in CM or CGM together with a low vertical resolution could 
result in "jumps" in the grids of max(F conv /F tot ) as function 
of Ted/logg. That is due to the fact that the convective zone 
expands in total size and retreats from the upper/mid photo- 
sphere when going from high-r e ff/low-g to low-r e ff/high-g. 
A change of the vertical extension of the convective zone even 
of a single layer can have a significant impact on the value of 
A(z), and hence on F com (~ A 8 in the low efficiency regime, see 
Gough & Weiss 1976 and references therein). Another aspect 
of the same problem was detected in computing FST stellar 
models with 72-layers atmosphere models as BCs. The thick- 
ness of the convection zone in the atmosphere A atm = z(F conv = 
0) - z(j = T p h) as a function of log g and was far from be- 
ing smooth, in particular in the low-r e if/low-g region of the 
HR diagram, where convection quickly extends into the higher 
atmosphere layers due to the contribution of H and H2 to the 
opacity. So, the resulting evolutionary paths in the HR diagram 
showed oscillations which disappear when new models with 
288 layers are used. 

The comparison between the 7\t)'s relations for CGM- 
atmosphere models integrated with 72 and with 288 layers (see 
Fig. 1 in Montalban et al. 2002) illustrates a result which from 
a different point of view was also discussed by Heiter et al. 
(2002a): in the case of ionization regions it is possible to cor- 
rectly follow the rapid increase of temperature only if the num- 
ber of layers is large enough. The same figure also points out 
a large uncertainty introduced for F and G stars (depending on 
their surface gravity), when the boundary condition is taken 
at optical depths r p h larger than 20, because of the onset of 
(nearly) adiabatic convection in the envelope of cool stars com- 
pared to the (almost) radiative temperature gradient of hotter 
stars. 

Atmosphere models are usually built to generate synthetic 
spectra to be compared with the observations. As a conse- 
quence, the number of mesh points per unit optical depth must 
be kept large in the region where the continuum and the lines 
are formed, that is at low optical depth, while few points are 
computed in between t = 10 and t = 100, though this is the 
recommended region (Morel et al. 1994) to adopt the values 
of temperature and pressure as BCs for the internal structure 
computation. The differences between T(t)'s relations corre- 
sponding to 72-CGM and 288-CGM atmosphere models show 
that low resolution models could induce non-negligible errors 
in the BCs, especially if r p h=100 is chosen as match point be- 
tween the atmosphere and the interior. 

3.3. Convection modeling in the interior 

For the interior as well, models are computed either by MLT 
or FST. If we want to build complete stellar models using the 
MLT treatment of convection and fit the Sun, we should use 
a-values in the interior (ai nt ) much larger than that one in the 
atmosphere (a a tm = 0.5) (see Sect.|4j. 

The MLT used in the interior strictly follows the formu- 
lation of Cox & Giuli (1968), who chose the V/A relation in 
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Fig. 3. FST evolutionary tracks for 1 M Q for two different 
choices of the tuning parameter a* and two different choices 
of the sub-photospheric boundary, T p h=10 and 100. The hori- 
zontal line refers to the models in Fig.|4] 

order to obtain numerical agreement with the results of Bohm- 
Vitense (1958). In the atmosphere, as mentioned above, the 
MLT formulation follows Castelli (1996). There, the choice of 
y and V/A parameters is such that for the same or-value, the ef- 
ficiencies of convection in Cox & Giuli formulation (Fcg) and 
in Castelli (Fc) hold the relation Fcg = 2/3Fc- This difference 
is much smaller than that introduced in the efficiency of con- 
vection by varying the a-value from -say- 1.7 to 0.5. 

The FST model of convection adopts the CGM fluxes, and 
the scale length at the radial coordinate r is the harmonic av- 
erage between z to p (distance to the top of convection increased 
by a*H plop ) and Zbot (distance to the bottom of convection in- 
creased by a*H p fr ot ), namely 

A = Ztop Zbat I faop + Zbot ) ( 1 ) 

This choice ensures that, close to the boundaries, z ~ Zi op or z ~ 
Zbot, and, far from boundaries, A ~ H p , as it should be expected 
(for details see Ventura et al. 1998a). For the models in which 
convection penetrates in the atmosphere, the quantity z to p has to 
include the distance from the bottom of the atmosphere (z(r p h)) 
to the top of convection in the atmosphere (z(ro)), and // Pjtop is 
also evaluated at z(ro) in the atmosphere. So, 

ZtopW = [z(t ) - z(Tph)] + CRmax - r) + a*H p (T ), (2) 

where r is the radius of the layer and R m ^ x is the radius of the 
last point computed by the integration from the interior, that is 
the radius of the boundary point (r p h). 

The fine tuning parameter a* for the computation of the 
interior can also be changed in order to fit the Sun. We will 
call ff* nt the parameter used for computation in the interior, in 



Fig. 4. Match between FST internal structure (a* =0.1, solid 
line; a* =0.2, long-dashed line) and the corresponding FST at- 
mospheres (short-dashed and dotted lines) for two different 
evolutionary phases of 1 M Q . The models correspond to the 
square points in Fig. [3] 

contrast with the fixed value a* atm = 0.09 adopted in the com- 
putation of the grid of atmosphere models. 5 

4. Effect of T ph on the evolutionary track location 

Since it is not reasonable by now to include a full atmosphere 
calculation into the computation of a whole stellar model, one 
separates the internal structure calculation, where the diffusion 
approximation is used, from the outer part. The calculation of 
the atmosphere is replaced by the values of the relevant phys- 
ical quantities at the optical depth r p h. This r p h must be large 
enough for the diffusion approximation for t > r p h to be valid, 
and small enough for the approximations adopted for the de- 
scription of the physics in the atmosphere to hold. Morel et al. 
(1994) proposed r p h > 10, depending on T e ff and logg of the 
model. 

If the calculations of the interior and the atmosphere were 
physically consistent, the stellar models computed with differ- 
ent choices of r p h should be equivalent. In practice, the physical 
inputs such as chemical composition, opacity tables, equation 
of state, and convection model used in the atmosphere are usu- 
ally not exactly the same as those used in the interior. 

The atmosphere models by Heiter et al. (2002a) adopt the 
opacity and the equation of state from Kurucz (1993, 1998), 
while our stellar evolution code considers the micro-physics 

5 We will approximate this value tt* lm = 0.09 with 0.1 in Table 
1 and in the text. This is due to the fact that we have used orL = 
0. 1 in the interior computations. For any purpose, 0.09 and 0. 1 provide 
practically the same results. 
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Table 1. Solar-composition computed models 
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interior 




ATLAS9 MLT 
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MLT 
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= 2.3 


10* 


ATLAS 9 MLT 
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= 0.5 


MLT 
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= 6.3 


100 (2) 


ATLAS 9 MLT 


a 
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MLT 
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1(4) 


ATLAS 9 MLT 
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MLT 
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= 1.85 
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ATLAS9 FST 
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FST 


a* 


= 0.1 


1,3,10, 100 


ATLAS9 FST 


a* 


= 0.1 


FST 


a* 


= 0.2 


1,3,10« 100 


AH97 MLT 


a 


= 1.0 


MLT 


a 


= 1.0 


3, 10, 100 


AH97 MLT 


a 


= 1.0 


MLT 


a 


= 1.9 


3, 10, 100 (5) 



FST cf=0.2, Heiter et al.02. t =10 



MLT a = 1.9, a ot =1.0 AH97. t =100 




3.75 3.7 3.65 

log T e „ 



Fig. 5. Several evolutionary tracks for 1 M Q , solar composi- 
tion, and MLT treatment of convection in the interior and in 
the atmosphere. Different values of a for each choice of r p h are 
required to fit the Sun with MLT atmospheres from Heiter et 
al. (2002a). Besides, a grey MLT model with the same opacity 
and EOS, and the new FST track with a* = 0.2 in the interior 
and FST atmospheres by Heiter et al. (2002a) are plotted. The 
thick line corresponds to a grey MLT model computed with the 
a value given by the 2D numerical simulations calibration of 
Ludwiget al. (1999). 

inputs described in Sect. 13.11 Fig. [5] shows, however, that, if 
the same convection model (CGM, a* = 0.1) is used in the 
atmosphere and in the interior, the HR diagram location of 
1 M evolutionary tracks is almost independent of the choice of 
T p h. Fig.[4]shows the structures (log T vs. log P) of the interiors 
and the corresponding atmospheres of 1 M G models at differ- 
ent evolutionary phases (PMS and Post-MS models marked in 
Fig-E]by square points). We see that, in spite of differences in 
the opacity tables and in the equations of state, the atmosphere 
structures between r = 10 and 100 almost overlap the values 
(log T vs. log P) obtained for the outer layers by integrating the 
internal structure. The differences introduced in T e ff due to the 
different choices T p h=10 or 100 are smaller than 0.5% along 
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Fig. 6. Stellar structure, interior and atmosphere for the so- 
lar model corresponding to the labeled evolutionary tracks as 
obtained with the evolutionary tracks plotted in Fig. [5] MLT 
models: atmosphere model with a- dlm = 0.5, r p h=100 (long- 
dashed line), T p h=10 (solid line), r p h=l (short-dashed-dotted 
line); atmosphere model from AH97, a atm = 1, r p h=100 (short- 
dashed line); FST model (long-dashed-dottedline). The circle, 
squares, triangle and cross indicate the match point between in- 
terior and atmosphere. 

the evolutionary track of 1 M Q with solar metallicity. This is 
not true if we consider smaller values of t^: for r p h= 1 or 
T p h= 3 the differences in temperature gradients between the in- 
terior and the atmosphere are larger, because the diffusion ap- 
proximation of radiative energy transport used in the internal 
computation is not valid. The temperature gradient computed 
in the stellar structure code begins to deviate from the atmo- 
spheric one at an optical depth lower than ~ 7, but the exact 
point depends on T e ff and on log g. Consequently, the details in 
the equation of state and opacities are of much less a concern 
for matching atmosphere and interior models than the region 
of validity of the diffusion approximation and the modeling of 
convection, to which we now turn our main attention. 

4. 1 . Solar calibrations 

It is generally accepted that any stellar modeling must be able 
to fit the Sun. Two parameters control the Sun's location in 
the HR diagram: the helium content which is determined by 
the solar luminosity, and the jump of specific entropy {As) 
between the photosphere and the adiabatic convection region 
(e.g., Christensen-Dalsgaard 1997) which determines the stel- 
lar radius. 

Ai= c p (V- V ad )Jlnp, (3) 

Jlnpo 
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where po refers to the pressure at a given point in the photo- 
sphere, and p ac j is the pressure at the point in the convective 
zone where (V - V a d) <k 1 . 

In MLT, (V - V a d) is directly related to the a parameter, 
and thus varying a allows the fit of solar radius. Hence, de- 
pending on the physical inputs, grey models provide values of 
a between 1.5 and 2.2 (Ventura et al. 1998a). Using the CGM 
model of convective flux and scale height, only a small varia- 
tion of (V - V a d) can be obtained by varying the tuning param- 
eter a* and for a solar grey model a* ^0.1 (CGM, Ventura 
et al. 1998a) is required. For the present models, using the re- 
lations T(t) provided by the grids of non-grey atmospheres, 
the value of a (and a* in the FST, but see later) required to 
compute the convective envelope fitting the Sun depends also 
on the choice of r p h. Eq. Q can be divided in two terms, 
Asatm(Tph) (jump of specific entropy from the surface to T p h) 
and Asi nt (r p h) (s jump between r p h and a layer deep enough 
where (V - V a d) ~ 0). Thus, the contribution of the model at- 
mosphere to the entropy jump As depends crucially on r p h and 
for any chosen r p h, As a tm(i" p h) must be smaller than As for a 
solar model to fit the Sun. 

We have considered several combinations of convection 
treatment in the interior and in the atmosphere, and several 
choices of r p h- Table 1 lists the different sets of models com- 
puted for the solar composition. The set of parameters fitting 
the Sun are indicated with a numerical super-script. The value 
of this index corresponds to the labels in Figs. [5] and |6] 

For CGM models, with the non-grey BCs the choice of 
a* = 0.1 as derived by the solar calibration with grey BCs 
provides a Sun ~ 85 K cooler than expected. To fit the Sun 
requires, hence, a larger a* value in the computation of the in- 
terior: viz. a* ~ 0.2 if r p h<10 (Fig.|4}. Models computed with 
consistent a* in the interior and in the atmosphere would have 
an intermediate value of a*. As shown in Fig.|4] an increase 
of a* by a factor 2 implies a lower temperature gradient in the 
layers deeper than r p h=10 and hence a higher T e g. The term 
a*H p (T()), however, affects only the scale length close to the 
boundary of the convective zone. Since the efficiency of FST 
convection rapidly increases downwards the convective zone, 
the atmosphere at r = 100 is quasi-adiabatic (see Fig. [2) and 
hence even a value of a* larger, but restricted to the interior, 
cannot increase the entropy in the adiabatic region. This ex- 
plains what is meant by saying that the FST model with the 
scale length as in Eq. or equivalently with the original 

suggestion from Canuto et al. (1996), is less parametric than 
standard MLT: if an FST-based solar model fails to fit the Sun, 
there is less possibility of achieving anyway a match by adjust- 
ing a", as this parameter only provides some fine tuning. 6 

Concerning the MLT calibration, Fig. [2] shows an atmo- 
sphere with T e ff and logg close to the solar values. The 
over-adiabaticity is quite high up to t = 100 and therefore 
As a tm(T"ph = 100) » Asatm(i"ph = 10). To keep As in the range 
required to fit the solar radius, the contribution of the term 
Asint(Tph) is controlled by increasing since (V— V a( j) ~ a -4 ^ 3 

6 In principle, the FST convective fluxes could also be used with a 
scale length aH p , but no grids of model atmospheres are available for 
this case. 



in the high efficiency limit. This explains the very different val- 
ues of aim required with different r p h's (a mt = 2.3 for r p h=10, 
and ai nt = 6.3 for r p h=100). Obviously, this results in a kink at 
the connection point, inducing discontinuities in the derivatives 
of the temperature gradient and in the over-adiabaticity (cf. 
Fig. |6j. A temperature structure similar to that obtained with 
FST could be obtained with MLT and a variable length scale, 
as in Schlattl et al. (1997). These authors found that between 
t(T = T e ff) and t = 20 a MLT 1-D atmosphere, computed with 
a - 0.5, agrees very well with the 2D-hydrodynamical model 
from Freytag et al. (1996). Hence they decided to use a model 
atmosphere computed with a = 0.5 down to r = 20 and intro- 
duced a mixing-length parameter a = a(j) whose dependence 
on t was chosen to model the temperature-pressure stratifica- 
tion of the 2D-model, and the parameters to fit the solar values 
of T e ff and luminosity (see their Fig. 1). 

Clearly, such procedures can be applied to the solar case, 
but how to extrapolate them to other regions of the HR dia- 
gram? Although these MLT-based models are built in such a 
way that the atmospheric convection prescription fits the solar 
Balmer lines better than the standard models with large a, and 
the internal structure integration fits the present solar radius, 
we have to warn about the poor physical meaning of extending 
this procedure to structures different from the Sun. To fit the 
Sun, the low a = 0.5 needed in the atmosphere is compensated 
by the larger value adopted in the interior, but the effect on the 
evolutionary track of a different a-value below r p h, is not the 
same over the whole HR diagram. As shown in Fig.|5] very dif- 
ferent PMS evolutionary tracks (by r e jf and shape) are obtained 
with parameter sets all fitting the HRD location of the Sun. We 
note that the PMS track of the 1 M model, obtained integrating 
the interior up to r p h=l, is, as expected, quite close to the MLT 
grey model. Moreover the coolest low-PMS tracks are obtained 
when the integration of the atmosphere (with low efficiency 
convection, a. dtm - 0.5, or 1.0) goes down to T p h=100. But, 
which are the good PMS tracks: those obtained with T p h=10, 
or with T p h=100? The convection properties of stellar layers at 
T p h are not the same all over the HR diagram. The fraction of 
the over-adiabatic region above r p h depends on T e g and log g. 
So, there is no physical reason in choosing an efficiency of con- 
vection below T p h and another above that point. The only jus- 
tification for doing so could be the match of a larger set of ob- 
served data, but this procedure would require calibration stan- 
dards throughout the HR diagram. On the other hand, the pro- 
cedure of mixing very different treatments of convection, based 
on solar calibration, finally results in physical inconsistencies, 
introduces a large dispersion (~ 200 K) in the low gravity re- 
gion of the HR diagram (see Hig-EJ, and reduces even more the 
predictive power of MLT based stellar modeling. 

Ludwig et al. (1999) proposed a calibration of the 
a(T e ff, logg) providing the same specific entropy jump in a 
grey atmosphere as their 2-D numerical simulations, and there- 
fore, the same stellar radius (although this does not imply the 
model will also have a realistic sub-photospheric stellar struc- 
ture). The thick line in Fig. [5] was computed with the same 
EOS and opacity tables, but adopting a grey atmosphere and 
MLT with the a(T e ff, \ogg) 2D-based calibration. We see that 
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for log L/L Q < 0.75 this track almost overlaps that obtained by 
non-grey FST modeling. 

Fig. [5] shows how different the sub-atmospheric layers of 
the "Sun's" can be, as obtained from five sets of parameters fit- 
ting the solar T e s and radius. Thanks to helioseismology it is 
possible to get information about the interior of the Sun, and 
several authors (Baturin & Miranova, 1995; Monteiro et al. 
1995; and Heiter et al. 2002b from non-grey MLT and FST 
models) have concluded that FST solar models yield sound 
velocities as a function of depth, which are in better agree- 
ment with the helioseismological data. Schlattl et al. (1997) 
as well, with their spatially varying mixing length parameter, 
obtained clear improvements with respect to previous compu- 
tations. It has to be noticed that even if the abrupt change of 
slope in the temperature profiles #2 and #3 of Fig. [6] is ade- 
quately smoothed, the frequencies resulting from both kinds of 
sub-atmospheric structures (FST and MLT with r p h=10) will 
be different and their differences can hence be probed by helio- 
seismology. 

The matching of model atmospheres on top of interior 
structure calculations raises the question if it is possible to put 
a constraint on the optical depth where convection becomes ef- 
ficient and if low values of a in an MLT treatment can be ex- 
cluded. The case of stars with shallow — in terms of H p — 
surface convection zones, such as A-type main sequence stars, 
can be described in terms of inefficient convection through- 
out the convective zone, from an observational point of view 
(Smalley & Kupka 1997, Heiter et al. 2002a, Smalley et al. 
2002), from the results of non-local convection models with- 
out a mixing length, and from numerical simulations (see the 
discussion in Kupka & Montgomery 2002, and Freytag 1995). 
Unfortunately, for stars with deep envelope convective zones 
the interpretation of the observational data remains ambigu- 
ous (cf. also Castelli et al. 1997, Heiter et al. 2002a) and the 
comparison with numerical simulations provides no simple al- 
ternative, i.e. a parameterization of temperature and pressure 
as a function of depth which could simultaneously fit both 
the stellar radii and the spectroscopic and photometric data. 
Obviously, the calibration problem becomes even worse for the 
case of evolutionary tracks running through extended regions 
of the HR diagram. In comparison with MLT, the FST based 
models offer the advantage of agreeing with more observational 
data without the need of drastically changing the scale length 
at some ill-defined matching depth, which, in the end, intro- 
duces at least another free parameter. Nevertheless, the discrep- 
ancies of these models with observed color photometry for the 
Sun (Smalley & Kupka 1997 and references therein) and with 
temperature gradients as obtained from various 3D large eddy 
simulations of solar granulation call for urgent further improve- 
ment in convection modeling. 

5. Convection treatment and PMS location 

From Fig.|3]it is clear that HR diagram location of the Hayashi 
tracks is strongly model-dependent, and the models with dif- 
ferent convection treatment — even selecting only those mod- 
els fitting the Sun — provide tracks of very different shapes and 
r e ff ranges (see also the first discussion given by D'Antona 



and Mazzitelli 1994). In this section we show the origin of 
differences and similarities among available PMS evolutionary 
tracks, often used in comparing and interpreting observational 
data. 



5.1. FST based models 

In this section we analyze the effect of different formulations 
of the FST model (different a*, plus the grey model) on the 
PMS evolution. As seen in the previous section, the solar cal- 
ibration with non-grey BCs requires a tuning of the parame- 
ter a* larger by a factor of two with respect to that used for 
the atmospheric grids. No feature from spectra and colors can 
constrain this parameter, but the deep structure of the atmo- 
sphere will be slightly different. The new value a* - 0.2 re- 
quired when r p h<10 is anyway within the limits given for stel- 
lar "overshooting" (Maeder & Meynet 1991). In principle, an 
inconsistency is introduced in this way between the treatment 
of convection in outer and internal layers. The effects of this in- 
consistency on the HR diagram are small, but different for each 
(T e ff, logg). For T e ff< 4700/f the atmospheric convection be- 
gins at increasingly shallower layers and becomes almost adia- 
batic deep in the FST atmosphere (see curve corresponding to 
?eff=4400 K and log g=3.8 in Fig.|2). This explains the small 
effect of different a* on the PMS tracks location in right side 
of the HRD compared with the effect on the PMS tracks corre- 
sponding to the highest masses considered. 

Fig.[7]shows the evolutionary tracks for masses from 0.6 to 
2.0 M obtained using: i) exactly the same treatment of con- 
vection in the atmosphere and inside the star (dotted lines), i.e. 
CGM fluxes and scale length with a* = 0.1; and ii) a slightly 
different scale length with ff* tm = 0.1, and a* = 0.2 below 
T p h=10. The effect on T e ff of changing a* by a factor 2 will 
depend on the quantity a = (A atm + Ai nt )/(a*H p ) (where A atm 
is the extent of the convection zone for t < r p h and A; nt is the 
distance to r p h). If a is much larger than one, the variation of 
r e ff will be negligible as the convective flux and the tempera- 
ture gradients will not change. On the contrary, if it is of the 
order of one or smaller, the effect on the temperature gradi- 
ents for the last points of the structure will be significant and 
as the gradients decrease, the r e ff's increase. The differences 
of T e ff, at constant luminosity, induced by the tuning of a* go 
from ~ 30 K at the lowest masses, to ~ 80 K for 1 .5 M , and 
~ 120 K for 2.0 M Q . This implies that, at given L and M, the 
maximum uncertainty introduced in T e ff is ~ 2% (< 1% for 
1 M Q ); and, at a given L and T e e the maximum uncertainty in 
mass is ~ 17%. 

Fig.|8]shows the effect of non-grey BCs on the location of 
the Hayashi tracks. DM97, 98 models were computed with the 
CGM formulation of convection in the interior and with grey 
BCs. The authors approached the problem of convection pen- 
etrating in the atmospheres of low-gravity/low-temperature 
models, integrating the interior until the top of the convective 
region. The present "complete" (and more physically coherent) 
FST models are much cooler than the grey ones, ~ 270 K for 
the lowest computed masses. 
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Fig. 7. FST evolutionary tracks with ATLAS9-FST atmo- 
spheres with a* =0.1 
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Fig. 8. Comparison between FST grey models (DM97), and 
"complete" FST models with the new ATLAS9 atmosphere 
models 

5.2. Effect of different MLT prescriptions 

We computed as well two grids of stellar models with the same 
micro-physics in the interior, but using NextGen atmosphere 
models (AH97) as BCs. NextGen includes a large number of 
molecular lines which, being the main source of opacity for 
low-temperature models, must have a significant role in the on- 
set of convection. These models use MLT with a = 1 .0 and pro- 



Fig. 9. MLT PMS tracks obtained using as BCs the atmosphere 
models by AH97 (a - 1), with two different values of a in the 
computation of the sub-atmospheric convection and with the 
match point as in BCAH98, T ph =100. 



vide super-adiabatic gradients smaller than those from the new 
ATLAS9 models. As a consequence, adiabaticity is reached at 
lower depths than in ATLAS9 (MLT a = 0.5) models and there 
is a smaller difference in choosing either r p h=10 or r p h=100. 

Anyway, the value = 1 is too small to fit the solar ra- 
dius, for which a mt = 1.9 is required. Fig. |5] shows the tracks 
computed with NextGen atmospheres (a atm = 1.0) down to 
r P h=100, and with a mX = 1 .0 (dash-dotted curve) and ai nt = 1.9 
(solid line) for all layers below. These curves are equivalent to 
those published by Baraffe et al. (1998, BCAH98). We point 
out that there is a large difference in T e ff (~380 K in MS and 
440 K in PMS for a 1M model) and that the shape of the 
Hayashi tracks as well changes with the a-value adopted in the 
interior. The NextGen tracks with the parameter fitting the Sun 
(ffatm = 1.0 down to r p h=100 and ai nt - 1.9) are cooler than 
ATLAS9 tracks with 

fatm — 0.5 down to T p h— 10 and a mi — 2.3. 
But Fig. |3 shows that T e s values as low as in NextGen mod- 
els can be obtained using the ATLAS9-tracks with a. dim - 0.5, 
T p h=100 and a\ nl = 6.3. The reason is just that, when we inte- 
grate the atmosphere down to r p h=100, we adopt a lower effi- 
ciency of convection in a larger region, and so the r e ff of the 
model decreases. 

We computed stellar models for AH97 BCs with a mt - 1 
and ffi nt = 1.9, for three different choices of r p h (3, 10 and 
100). As in the FST case (Sect.|4}, when we use a mt - 1.0, the 
same as in the atmosphere models, there is only a small differ- 
ence between evolutionary tracks whose interior has been in- 
tegrated until r p h=3, 10, or 100. The difference between r p h=3 
and r p h=100 tracks, due to the incoherence introduced by the 
diffusion approximation used in the interior integration, is of 
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Fig. 10. Comparison between MLT models using ATLAS9 at- 
mospheres and AH97 atmospheres until r p h=3. In both cases 
the parameters are chosen to fit the Sun. The dashed-line is a 
1 M evolutionary track computed with a m = 1.9 and with 
AH97 atmospheres until r p h=100. 



the order of 40 K in the low-gravity domain. When we use 
QTint = 1.9, but the atmospheric grids still have a — 1.0, the 
PMS track location depends as well on r p h. Thus, the r p h=3 
tracks are now ~ 170 K hotter than the T p h=100 tracks (see 
Fig.^|for the relative location of 1M tracks, thick-solid line 
and dashed line). 

In Fig.^|we plot the evolutionary tracks calibrated on the 
Sun and integrating the internal structure up to r p h=3 for both 
atmosphere models, NextGen and new MLT ATLAS9 mod- 
els. The overlap of both sets of models for the whole common 
T e ffflogg domain is, at first sight, surprising: it points out the 
similarity of the opacity in the atmosphere models (NextGen 
versus ATLAS9) in this domain. Generally, at t = 3 the gra- 
dient of temperature still follows the radiative gradient of tem- 
perature, and hence the convective fluxes and the choice of a atm 
do not have a relevant role. 

These results confirm that for r e ff>4000 K the convection 
model is much more important than the contribution of addi- 
tional molecules to opacity already included in ATLAS9 (as 
concluded by DM94 on the basis of grey models). 

5.3. FST versus MLT 

The main difference between FST and MLT that we must keep 
in mind is that in the region of low efficiency FST is much 
less efficient, while, in the high efficiency region, the FST con- 
vective fluxes are ~10 larger than the MLT ones, so FST de- 
scribes a much more efficient convection in deep layers. This 
means that, deep in the atmosphere, the MLT models yield 
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log T e „ 

Fig. 11. PMS tracks obtained using parameters that fit the Sun, 
for FST and MLT treatment of convection using the new MLT 
and FST ATLAS9 atmosphere models 

higher temperature gradients than the FST ones (cf. Fig.|6]for 
log(P) > 5.3), and hence lower r eff . DM97, 98 showed that, 
when grey BCs are used, the difference in T e g between FST 
and MLT tracks for solar mass tracks fitting the Sun was in be- 
tween 5-10% (~ 250 K for 1 M Q ). This difference is relatively 
small in absolute value, but quite large, if we consider that PMS 
tracks of different masses are nearly parallel and rather close to 
each other. Thus, the masses of TTauri determined by compar- 
ison between their location on the HRD and either the MLT or 
FST evolutionary tracks, could differ up to 50%. 

If non-grey BCs are used, the T e s location of MLT tracks 
also depends on r p h and therefore a direct comparison with 
FST tracks is not obvious. Fig. shows the new FST and 
MLT Hayashi tracks computed with the parameters {a* and 
a) calibrated on the Sun and T p h=10. We note that the differ- 
ences increase as the stellar mass decreases, and they are almost 
negligible for M > 1 M . There are two reasons that explain 
why the MLT and FST tracks are now so close to each other: i) 
Fig- D] shows that there is a domain in (T e g, \ogg) where both 
sets of atmosphere models with low convection efficiency pro- 
vide very similar atmospheric structures for t < 10, and hence 
similar BCs. ii) The low efficiency of MLT convection below 
t = 10 is more or less compensated by a higher value of a 
(2.3) in the interior. The total effect is such that the largest dif- 
ferences in r e jf are of the order of 100 K (~ 2%), that induce 
an uncertainty in the mass determination of observed PMS stars 
between 10-17%. 

Note also the difference between FST and MLT MS models 
at ~ 6500 K, where convection in the atmosphere is highly in- 
efficient (this can be expected from their different atmospheric 
temperature structure as shown by Fig.^ see a ls° D' Antona et 
al. 2002). 
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Fig. 12. PMS tracks of 2.0, 1.0, and 0.6 M stars for three dif- 
ferent chemical compositions ([M/H] = -0.3, 0.0, +0.3). All 
these tracks have been computed using as boundary conditions 
the new ATLAS9 (CGM, a* = 0.1) atmosphere models at 
T p h=10 and the CGM formalism for the convection in the in- 
terior with a* = 0.2. 

6. Metallicity and PMS location 

We have computed FST-based models also for two other metal- 
licities, [M/H]=-0.3 and +0.3. All of them included the CGM 
convection model with a* =0.1 until r p h=10, and a* = 0.2 for 
the layers below. Fig. 1121 shows that changes of the chemical 
composition by a factor 2 with respect to the solar one pro- 
duce a small effect on the T e g location of PMS tracks for FST 
models: ~ 100 K hotter for [M/H]=-0.3 and ~ 100 K cooler 
for [M/H]=+0.3. In Fig.^]we compare our FST models with 
MLT models by Siess et al. (2000) for 1 M . At log L/L G = 0.6, 
the FST-based model with [M/H] =-0.3 is ~ 130 K hotter than 
that one with solar composition, while MLT models differ by 
~ 170 K. Note that in the Siess et al. (2000) models the Z=0.02 
solar track does not fit the solar T e g. Nevertheless, the Siess 
et al. solar metallicity (Z=0.02) PMS models, being cooler, do 
not make a bad job in reproducing the observations, according 
to the analysis by Simon et al. (2000). 

7. Comparison with the observations 

The theoretical uncertainties in the location of the PMS, con- 
firmed in the present computation, are such that direct mass 
determinations of PMS are needed to constrain the models. 
In particular, the eclipsing double-lined spectroscopic binaries 
allow simultaneous determination of masses and radii of the 
components. Only one binary of this type is known by now 
to contain two PMS components, namely RXJ 0529.4+0041 
(Covino et al. 2000) while several other PMS spectroscopic 
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Fig. 13. Evolutionary tracks for 1 M Q for metallicity Z = 0.01 
and Z = 0.02 (Z ) and two treatments of convection: FST with 
a* = 0.2 and MLT from Siess et al. (2000). 

double-lined non-eclipsing binaries (Melo et al. 2001) can just 
constrain the minimum mass of the components. In Fig. [2] 
we plot the masses derived by comparing the observationally 
derived T e ff (and luminosity) with the theoretical evolutionary 
FST-tracks versus the dynamical masses. The error bar on the 
theoretical mass reflects the observational error on the temper- 
atures quoted by the authors (~ 100 - 200 K). A very small 
error is generally attributed to the dynamical masses. Apart 
from the case of RXJ 0529.4+0041, for which the inclination 
i is known, the plotted dynamical mass is the minimum mass 
Msin 3 i. Many points remain below the diagonal line in Fig. 1141 
and the discrepancy could be larger if we had a precise indi- 
cation of an i < 90. This result, although still preliminary in 
view of the very few observational points and of the uncertainty 
of T e s determination of the PMS objects, indicates that even 
with the inclusion of non-grey FST model atmospheres, the 
FST convection is too efficient to describe the PMS. Moreover 
it confirms the preliminary comparisons based on grey FST 
models by Simon et al. (2000) and Covino et al. (2001). It is 
possible that the error bars on T e ff as provided by the authors 
were underestimated. In order to make the mass derived for the 
coolest component of RXJ 0529.4+0041 compatible with ob- 
servations, the error bar should grow up to ~ 400 K. Since 
the PMS-FST tracks are very close to each other in T e $ (at 
log L/L Q = 0.3, Ar eff between 1 .5 and 0.6 M Q is only -650K), 
such a large error in T e g also for the other stars would imply, 
however, that no useful constraint on the models could be ob- 
tained from this comparison. 

Baraffe et al. (2002) found that most masses in PMS bi- 
naries are better reproduced with = 1 .9, but some others 
require a mt = 1, while in principle there is no physical reason 
to justify this difference in the efficiency of convection in the 
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Fig. 14. The mass determined from the HR diagram position 
with respect to the evolutionary FST tracks is plotted with 
respect to the dynamical mass for the double lined spectro- 
scopic and eclipsing binary RXJ0529. 4+0041 (full circles); the 
other points are for the double lined spectroscopic binaries 
RXJ0541. 4-0324 (empty squares); RXJ0532. 1-0732 (empty 
triangles); RXJ0530.7-0.434 (empty circle, only one because 
M,T e s and L are the same for both stars), in which the mini- 
mum dynamical mass is plotted. The data are from Covino et al. 
2001. In addition, the location of TYCrA (full triangle) (Casey 
et al. 1998) and GMAur (full square) (Dutrey et al. 1998) is 
shown. 



same region of the HR diagram. In fact, we have seen in Sect. 5 
(Fig-EJ that for a given mass, the MLT tracks may differ even 
by 500 K if we consider both models with small and large a and 
disregard any solar constraint. Furthermore, for solar calibrated 
models, a variation of r p h from 3 to 100 could produce a dis- 
placement of their PMS tracks of the order of 200 K (compare 
Fig.[9]with Fig. I10i. No progress in understanding is possible 
with this philosophy. In particular, there may be hidden param- 
eters which affect the tracks location that could not be isolated 
unless we can — at least in part — falsify the models. 

FST, contrarily, does not allow large displacements on the 
HRD by changing the parameters. The FST convection has 
also been shown to be very useful both to describe the solar 
structure (see the introduction) and many other different evolu- 
tionary phases (e.g. the DB white dwarfs pulsation boundary - 
Benvenuto & Althaus 1997-, the Hot Bottom Burning phase in 
the intermediate mass Asymptotic Giant Branch stars -Ventura 
et al. 2001 and 2002-, the fast transition between deep and 
atmospheric convection in main sequence stars -D'Antona et 
al. 2002): this may lead us to suspect that there is some other 
physical input specific of PMS models, which modifies convec- 
tion in the PMS stars. Just as an example, we notice that GM 
Aur, which has a dynamical mass estimation of 0.84 + 0.05 M Q 



(Dutrey et al. 1998) and an evolutionary mass of 0.5 + 0.2 M 
also shows a magnetic field of 4 - 5 10 3 Gauss (Johns-Rrull 
et al. 1999). D'Antona et al. (2000) have shown that the purely 
thermal effect of a magnetic field deeply modifies the behavior 
of atmospheric convection, precisely in the required direction 
of enlarging the temperature gradients and lowering the models 

8. Summary and Conclusions 

The new ATLAS9 atmosphere models based on the CGM treat- 
ment of convection have allowed us to compute, for T e ff > 
4000 K and three different metallicities, new grids of non-grey 
FST-PMS evolutionary tracks with the same convection treat- 
ment in the interior and in the atmosphere. They provide a sig- 
nificant improvement with respect to DM97,98 grey FST-PMS 
tracks. Furthermore, we have examined one by one the param- 
eters affecting the T e g location of PMS tracks analyzing the 
"traps" hidden in the modeling of convection in PMS and tak- 
ing advantage of the use of new grids of stellar atmospheres 
constructed with different assumptions for convection. 

Concerning the choice of the match point r p h between the 
atmosphere and the interior, we have pointed out that the differ- 
ences between the micro-physics (opacity tables and equation 
of state) inputs in the interior and the atmosphere models in- 
troduce a negligible uncertainty on the derived HRD location. 
Actually, the spatial resolution of the deepest layers of the at- 
mosphere model is much more relevant for the choice of T p h. 
We have seen that low resolution models could induce non- 
negligible errors in the BCs, especially if r p h=100 is taken as 
the match point between atmosphere and interior. This choice 
may otherwise appear more suitable due to the flatter temper- 
ature gradient and the ascertained validity of the diffusion ap- 
proximation of radiative transfer at this depth. 

We show that, while a successful solar calibration of a 1M 
model adopting a low efficiency of convection in the atmo- 
sphere can be achieved by adequately increasing the efficiency 
of convection in the interior, sets of parameters which allow 
obtaining a good solar radius do not provide the same PMS 
tracks. In fact, different couples of parameters (ai nt ,T p h) intro- 
duce a dispersion in the T e s of the PMS track of 1M Q of the 
order of 200 K (if we take into account also the FST track, the 
uncertainty on T e s grow up to 250 K). Since there is no physical 
reason to have a different efficiency of convection at each side 
of a fixed r p h, the introduction of another parameter reduces 
significantly the predictive power of such models: in non-grey 
stellar models computed using different convection in the inte- 
rior and in the atmosphere, the parameters are at least three: 
ami, ffatm and r ph . 

While in the MLT models the value a- dtm = 0.5 is chosen to 
fit the Balmer lines, in the CGM ones, the fine tuning param- 
eter a* cannot be constrained by the spectral features of solar 
type stars, since a* affects mainly the gradients in the deep lay- 
ers of the solar atmosphere. The value adopted from the grey 
solar calibration by Canuto et al. (1996) turns out to be too low 
for the non-grey model of the Sun, and a variation of a? . of 
a factor of two with respect to <?* tm is required. A solar model 
with exactly the same convection treatment in the interior and 
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in the atmosphere should adopt an intermediate value of a*. 
Nevertheless, by changing a* ml within this range we are not able 
to produce significant displacements over the HR diagram, and 
the Ar e ff on the PMS tracks due to this change is smaller than 
100 K for M < 1.5M . 

The evolutionary track using grey atmospheres and a MLT 
treatment of convection with a a(T e e, \ogg) value provided by 
the 2D-model calibration (Ludwig et al. 1999) is quite close 
(in a large log g /T e ff domain) to that one obtained with a non- 
grey FST model fitting the Sun. This implies that the specific 
entropy jump between the photosphere and the adiabatic zone 
is almost the same in the FST based models and in the 2D- 
hydrodynamical atmosphere models. However, one should ex- 
pect the specific entropy jump and hence PMS evolutionary 
tracks obtained from a 3D-hydrodynamical model atmosphere 
calibration to be different. As shown by Asplund et al. (2000), 
for solar surface convection numerical simulations in 2D yield 
higher temperatures in the deeper layers and lower ones in op- 
tically thin regions when compared to their 3D counterparts 
(see their Fig. 9). The calibration by Trampedach et al. (1999) 
which is based on such 3D numerical simulations cannot be 
used safely to calculate PMS evolutionary tracks, because the 
model calculations used to derive the coefficients of their fitting 
functions are limited to essentially the main sequence band. 
Ludwig et al. (1999) explicitly warn their readers that their fits 
quickly loose their meaning outside the region covered by their 
own model grid. Fortunately, their grids have been computed 
for a much wider region of the HRD which in fact is just large 
enough to include the solar PMS track as shown in Fig.[5]within 
the domain valid for their calibration. 

We obtain quite similar PMS tracks in some regions of the 
HRD by using FST or MLT (a atm = 0.5, a mt = 2.3, T ph =10) 
stellar models. The reason is that MLT (a = 0.5) and FST 
atmospheres have {P, T) structures which overlap for optical 
depths smaller than r ~ 10, and the higher efficiency of FST 
at larger optical depths is compensated by a larger value of the 
mixing length parameter. The maximum Ar e ff is of the order of 
100 K, for the lowest considered masses. 

The comparison between AH97 and the new ATLAS9- 
based models shows that: i) both sets of models are equivalent 
if the match point between atmosphere and internal structure 
is chosen not too deep in the atmosphere (e.g. r p h=3, Fig. I10i. 
Since in this region the temperature gradient follows the radia- 
tive gradient, and the radiative gradient is determined mainly 
by the opacity, the overlapping between both sets of models 
all over the T e ff -log g domain considered allows us to conclude 
that the improvement of opacity in AH97 models has no rel- 
evant role for the location of stellar evolution tracks in this 
region of the HR diagram; ii) the differences between MLT- 
ATLAS9 models and AH97 models are only due to the spe- 
cific treatment of the over-adiabatic convection. Therefore, at 
r e ff>4000 K the convection model is much more important 
than the contribution of additional molecules to the opacity 
already included in ATLAS9 (as concluded by DM94 on the 
basis of grey models). 

Examining some PMS binaries, which provide independent 
determination of the stellar masses, the new FST r p h=10 mod- 
els provide an on average too efficient convection with re- 



spect to the masses and HR diagram location of these binaries. 
Covino et al. (2001) and Steffen et al. (2001) found that the 
best fit is obtained when MLT with a = 1 (in the interior and 
in the atmosphere) is adopted. Since MLT based models must 
adopt a value of a larger than 1 to fit the Sun, and since FST 
treatment (that provides a quite efficient convection) is able to 
fit the solar radius, the comparison with observation suggests 
that convection in the PMS evolutionary phase seems to be 
less efficient than for the Sun. The FST model cannot be tuned 
too much, but it nevertheless provides very valuable results for 
many different evolutionary phases. Therefore, instead of dis- 
missing the model at this point, we suggest that there may still 
be another physical parameter affecting the location of PMS 
tracks in the HR diagram. Processes such as, e.g., interaction 
with a circumstellar disk (Flaccomio et al. 2003), acceleration 
and/or braking of rotational velocity during PMS evolutionary 
phase, or presence of a dynamo magnetic field (Ventura et al. 
1998b, D' Antona et al. 2000) could modify the convective tem- 
perature gradients in the outer layers of PMS stars. 

Finally, also the sub-photospheric structures of the Sun, as 
obtained by different choices of BCs and convection model- 
ing, are quite different. The helioseismological data allow now 
to discriminate between different solar structures and indicate 
that the Sun is better reproduced with the FST-based model 
than with the MLT one. However, there are other aspects of 
surface convection in the Sun that cannot be fitted with a lo- 
cal ID model such as FST, as discussed in Sect. l3.2.T1 and l4~T1 
and in references cited therein. Major progress towards non- 
local and non-homogeneous models of convection suitable for 
evolutionary track computations is thus urgently needed to im- 
prove current stellar models. 
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